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Abstract 

We investigate constraints on power spectra of the primordial curvature and ten- 
sor perturbations with priors based on single-field slow-roll inflation models. The 
Hubble slow-roll parameters are included in cosmological parameters and the primor- 
dial power spectra are generated using the inflationary flow equations. Using data 
from recent observations of CMB and several measurements of geometrical distances 
in the late Universe, we perform Bayesian parameter estimation and model selection 
for models that have separate priors on the slow-roll parameters. The same analysis 
is also performed adopting the standard parameterization of the primordial power 
spectra. We confirmed that the scale-invariant Harrison-Zel'dovich spectrum is dis- 
favored with more significance than previous studies. While current observations 
appear to be optimally modeled with some simple models of single-field slow-roll 
inflation, data is not enough constraining to distinguish these models. 



1 Introduction 



Cosmological observations, including the cosmic microwave background (CMB), large scale 
structure, baryon acoustic oscillation (BAO), type la supernovae (SN) offer opportunities 
to probe new physics far beyond the reach of experiments in terrestrial laboratories. One 
of such physics is inflation, that solves various problems in the hot universe scenario 
P~H5]. In addition, inflation also explains generation of initial perturbations for structure 
formation [6HT3]. At present, Inflation is an essential part of our best description of the 
Universe. 

The simplest class of models of inflation is so-called single- field slow-roll inflation [HE], 
where potential energy of a single scalar field (inflaton), whose field value varies slowly, 
drives an exponential expansion of the Universe. Gaussian, nearly scale-invariant primor- 
dial curvature perturbation can be generated from the vacuum fluctuation of inflaton in 
a quasi-de Sitter spacetime. Such primordial perturbation gives excellent fits to various 
data from cosmological observations, which makes single-field slow-roll inflation highly 
attractive. However we yet know quite little about inflaton and its potential. 

Identification of inflaton is of particular interest both in cosmology and particle physics. 
As a model of single-field slow-roll inflation left its vestige in the late Universe in an ob- 
servable way through generation of the primordial curvature and tensor perturbations, 
constraints on models have been mainly investigated through constraining power spectra 
of these primordial perturbations, V^(k) and Vh(k), using various cosmological observa- 



tions^]. One of the most familiar ways may be to adopt the standard parameterization 
for the power spectra, A s , r, n s , n t , a s , etc. (See Eqs. (ITMT9]) ). and derive constraints 
on these parameters. Some more involved analyses have also been performed focusing 
on reconstruction of the potential of the inflaton or slow-roll flow parameters [TH432] . 
In either way, deviation from the Harrison-Zel'dovich (HZ) spectrum (i.e. r = and 
n s = a s = ... = 0) allows us to probe models of single-field slow-roll inflation. Cosmological 
data now signifies some spectral features in the power spectrum of the curvature perturba- 
tion. When we assume a power-law for curvature perturbation spectrum V^(k) oc A s /c ns_1 
and absence of the tensor perturbation Vh{k) = 0, recent CMB data from WMAP give 
—0.065 < n s — 1 < —0.009 (95% C.L.), suggesting a significant deviation from the HZ 
power spectrum [33] . 

However, constraints on the primordial power spectra are highly dependent on pa- 
rameter spaces which we investigate. For instance, if running of the spectral index a s is 
included, the same data allow relatively large a s [33] ■ This kind of issues commonly arise 
when we try to constrain the primordial power spectra in parametric ways. Consequently, 
shapes of reconstructed power spectra also differ depending on the choice of parameter 
spaces, which finally affects selection among inflation models. 

Since we do not in advance know a parameter space where we should explore a con- 
straint on the power spectrum, we also need to examine whether the parameter space 



There are also several other probes for inflation models, including primordial isocurvature perturba- 
tions and non-Gaussianity in primordial perturbations. 
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is appropriate. More generally, appropriateness of a model, which possesses its own 
prior assumption, should also be discussed. A guiding principle in looking for an optimal 
model is Occam's razor, which penalizes unnecessary assumption in describing observa- 
tions. Bayesian model selection is Bayesian implementation of Occam's razor, which is 
now frequently applied in the context of cosmology [33H52] . In particular, the authors 
in Ref. [25] have adopted Bayesian model selection to assess optimal orders up to which 
reconstruction of inflaton potential should be performed. Moreover, in Ref. [53] Bayesian 
model selection is directly applied to distinguish some class of inflation models. 

Motivated by [25|,|53] and other earlier studies, in this paper we investigate an optimal 
constraint on primordial perturbation spectra and make comparison of single-field slow- 
roll inflation models using recent cosmological observations, which would be a subject of 
great interest [51]. For these purpose we make vigorous use of Bayesian model selection. 
We compare Bayes evidences for several models which have separate priors on inflationary 
slow-roll parameters, or parameters of primordial power spectra. Each of these models 
can be regarded as representing some class of single-field slow-roll inflation models. In this 
paper we implicitly assume a Freedman-Robertson- Walker universe and adopt natural 
units h = c = Mpi = 1. 

This paper is organized as follows: In Section [2] we briefly review the Hubble slow- roll 
flow equations, which are adopted in our analysis. Then in Section El some essences of 
Bayesian model selection are also reviewed. The main part of our paper is Section HJ where 
we investigate constraints on the primordial power spectra using data from recent obser- 
vations of CMB and geometrical distances in the late Universe. By employing Bayesian 
model selection, an optimal constraint and comparison of models of single-field slow-roll 
inflation are investigated. The final section is devoted to summary and future outlook. 

2 Hubble slow-roll flow equations 

We make use of the Hubble slow-roll (HSR) flow equations [16], [19], [55] , a specific version 
of more general inflationary flow equations jSHJEI]- Our algorithm is basically identical 
to slow-roll reconstruction presented in [SSIED] and we here quote some key consequences. 
In this paper, we assume models of single-field slow-roll inflation have canonical kinetic 
terms and do not investigate non-canonical models. 

During an epoch of slow-roll inflation, the inflaton field can be safely assumed as 
a monotonic function of the time and hence regarded as a generalized time coordinate. 
Instead of solving the Hamilton- Jacobi equation [57] for given inflaton potential, we rather 
solve a set of the HSR flow equations given by 




(for £ > 1). 




(2) 
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Here and hereafter we denote derivatives respective to as prime (e.g. H' = dH/dip). 
Notice that <fi plays a role of time variable in this formalism. These flow equations are 
solved once we specify the initial values for e# and £ Xh at a fiducial 0*, 



Ah* 



e X H ( 



(for £>1) 



(3) 
(4) 



If we choose fiducial HSR parameters £ \h* — for £ > M, the HSR parameters £ A#(0) 
for i > M also vanish at any and hence the flow equations are truncated at order M. In 
this case, the Hubble expansion rate can be exactly solved [EE]- Without loss of generality, 
we can choose 0* = and being a decreasing function of time (i.e. H'((p) > 0). Then 
we obtain 



H{(j>) = #*(! + B 1( j) + B 2 <j) 2 +-B 



M+l<- 



(5) 



where = H(<f>*) = H(0) and the coefficients of the Taylor expansion, Bi (i = 1, . . . , M + 
1), are given by the HSR parameters at the fiducial point, 



B 1 = a/ Ait e H *, 

(47T)' 



Be+i 



(£ + l)\Bf 



~ Xfj* (for£>l). 



(6) 
(7) 



Since we have specified the function H(<f)), we already know the dynamics of the back- 
ground universe and related variables as functions of <fi. For example, e-folding number 
N((f>) at time 4> is given by integrating 



dN 



Air 

elm' 



Similarly, a wave number that exits the horizon at time is given by integrating 



d In k 



Air 



e H (<f>)). 



(9) 



Once we fix a fiducial wave number fc* = fc(0#), there is a one-to-one correspondence in 
k and <fi by Eq. (Q, as long as slow roll of the inflaton does not break down. Then any 
other function of is rewritten as that of k. For later convenience, we choose the e- 
folding number N(k) so that it vanishes at a wave number k = 10~ 4 Mpc~ 1 , which roughly 
corresponds to the scale of current horizon ~ 10 Gpc. 

In this paper, we truncate the HSR flow equations at M — 2. This is because higher 
order HSR parameters may not be important due to the limited range of observable wave 
numbers O(10~ 4 ) < k < 0(0.1) Mpc" 1 and insufficient accuracy of data at present or in 
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the near future. The fiducial HSR parameters up to M = 2 are nothing but the usual 
slow-roll parameters e(0), 7/(0) and £(0) and they are given by 



<4>) 



^(0), 



(10) 

(11) 

(12) 



One of the most important prediction of single-field inflation models is that the pri- 
mordial curvature and tensor perturbations are not independent of but related to each 
other since both are generated from the dynamics of the single scalar field 0. Thus there 
exists the consistency relation for single-field inflation models. At second order in slow-roll 
approximation, the power spectrum of the primordial curvature and tensor perturbations, 
V^(k) and Vh(k), respectively, are given by [59] 



VAk) 



(2C+l)e(0) + CrK0)F 



7T6 



(0) 



0=0(fc) 



= i6[i-(c + i W )F , 



7T 



(13) 
(14) 



$=4>(k) 



where C = -2 + In 2 + 7 ^ -0.729637 and 7 is the Euler constant. Although Eqs. (I13til4p 
are originally derived assuming e and rj are constants [59], evolution of these variables are 
practically very small at observable scales k, which validates the use of Eqs. ( IT3TIT4"|l l#^l. 

The standard parameters of the primordial power spectra are given by Taylor expand- 
ing the logarithm of Eqs. 



n. = 



dink 

P h 



1 = 1 + 2rj* - 4e* - 2(1 + C)el 



3-5C 



3-C 



k=k, 



Pi 



16e*(l + 2C(e* -77*)), 



(15) 
(16) 

(17) 



d 2 lnV t 



O.S 



n t = 



din k 2 
d\nV h 



dink 



k=k t 



k=K 



-- -2& - 8e^ + 10e,77„ 
-2e*-(3 + C)e^ + (l + C)e,r/, 



(19) 



where C = 4(ln2 + 7) — 5. With these standard parameters, we can approximate the power 



# 2 Several authors [SSIEQ], have discussed effects of differences in the power spectra calculated from 
exact solution of the wave equation, approximation with slow-roll parameters up to the 2nd order (i.e. 
Eqs. (| 13til4[) ^1 . They conclude that differences among the exact and these approximated power spectra are 
not significant with a certain prior on the e-folding number N. 
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spectra by 



P c {k) = A s exp 
V h {k) = rA s exp 



(n.-l)ln— + -a s In- 



n t m — 



(20) 
(21) 



Along with the power spectra in Eqs. (I13til4p . the duration of the inflation is also 
important. We impose a prior on the e-folding number, iV > 25. This prior roughly 
means that the energy scale of inflation should be higher than TeV |23p^l . Imposition of 
this prior ensures that inflaton does not affect physics below electroweak scales. However, 
this does not necessarily mean that a given model of single-field slow-roll inflation with 
certain (e*, 77*, £*) is excluded if it predicts N < 25, for slow- roll parameters at higher 
orders can maintain slow-rolling of the inflaton further. Thus we cannot strictly restrict 
models of single-field slow-roll inflation with a prior iV > 25. In Section [4j as default, we 
adopt a prior N > 25 in investigating constraints on the primordial power spectra with 
inflationary priors. However, we also investigate constraints without the prior N > 25, 
which would be informative in discussing what are plausible models supported from data. 



3 Bayesian model selection 

As we have mentioned in Introduction, we adopt Bayesian model selection in order to 
compare different models of single-field slow-roll inflation. Before presenting details of our 
analysis, let us briefly review some essences of Bayesian model selection. We also refer 
to |l£lll71ll9] for more detailed reviews. 

Given data, a joint posterior distribution for a model M and its model parameter 0, 
P(0, M|data), is given by the hierarchical Bayes theorem, 

= P(da t a|e,M) f( e|M)PW (22) 
P(data) 

where P(data|G,M) = £(G) is the likelihood function, P(G|M) is the prior distribution 
for Q specified by M, and P(M) is the prior distribution for M. The remaining P(data) is 
an irrelevant normalization constant. Then the posterior distribution for M, P(M|data) 
is given by marginalizing P(0, M |data) over the model parameters 0, 



P(M|data) = j d0P(0, M|data) 
P(M) 



P(data) 



J d0P(data|0,M)P(0|M). (23) 



# 3 The role of a prior on the e-folding number is detailedly discussed in [60) 
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The final integral in Eq. ( )23|) is called Bayes evidence E(M), 

E(M) = J c/6P(data|6, M)P(6|M), (24) 

which measures marginalized likelihood of the model M. Thus the relative likelihood of 
different models M; and Mj is assessed by the ratio of their Bayes evidences, 

* = In ^ (25) 
13 E(Mj) v ; 

which is called Bayes factor. Bayesian model selection is implemented by estimating Bayes 
factor which measures relative likelihood between two different models. Jeffreys' scale is 
often adopted to connect numbers and semantics, which states: for \B{j\ < 1 the evidence 
is not significant; 1 < \Bij\ < 2.5 significant; 2.5 < \Bij\ < 5 strong; and 5 < \Bij\ 



decisive #4 



In usual parameter estimation, a posterior distribution of parameters for a model M, 
P(0|data, M) is investigated. However, our purpose in this paper of comparing different 
models of inflation is not achieved via P(©|data, M), where a model M is by definition 
assumed to be true. This is another way of representing our statement in Introduction, that 
different models cannot be compared with a fixed parameter space. Instead, comparison of 
different models is achieved by the Bayes factors B^. This is why we adopt the Bayesian 
model selection in our analysis. By assigning a prior P(Q\M) to reflect prediction of a 
class of single-field inflation models M, the likelihood of M is measured with the Bayes 
evidence E(M) and two different classes can be compared. 

On the other hand, we can also explore an optimal constraint on parameters by 
adopting Bayesian model selection. By comparing different models, which have sepa- 
rate priors on O, an optimal constraints can be obtained from a posterior distribution 
P(0|data, M) of the parameters for a model which has higher Bayes evidence than other 
models. 

One remark should be mentioned. The Bayes evidence explicitly depends on a prior 
probability P(Q\M) due to its normalization. For example let us consider that some two 
different top-hat priors are imposed on a same parameter. Even though the parameter 
region of high-likelihood is sufficiently covered by both of the priors, the resultant Bayes 
evidences differs, that are inversely proportional to the ranges of the top-hat priors. This 
is very contrastive to the usual parameter estimation, where the posterior probabilities 
from these two different priors do not differ much. 



4 Analysis with observational data 

Here we investigate constraints on the primordial power spectra and make comparison of 
single-field slow-roll inflation models. First of all, since there are quite a large number 

# 4 The choice of numbers and semantics is more or less arbitrary and differs literature by literature. 
For instance, a more conservative statement is suggested in [46] . 



6 



of theoretical models of single-field slow-roll inflation, it is difficult to analyze them indi- 
vidually. Instead we consider several different models which have separate priors on the 
primordial power spectra. Each of these models can be regarded as representing some 
class of single-field slow-roll inflation models. We adopt two parameterizations of the pri- 
moridial power spectra. One is the standard parameterization of Eqs. (120112 ip . that has 
A a , n s , r and a s . As we are interested in single- field slow- roll inflation models, we impose 
the standard inflation consistency relation, n t = — r/8, when we adopt the standard pa- 
rameterization. The other is the HSR parameterization that has A s , e*, 77* and £*. By 
adopting these two parametrizations, the models we compare are listed in Table [TJ The 
reference model is Mrz that has the HZ primordial power spectrum. M ns , M nsr and M ns0ls 
have different top-hat priors on the standard parameters n s , r and a s , that are described in 
the last column of Table [H On the other hand, M e , M v , have different top-hat priors 



on the HSR parameters e* and 77* that are also described in the last column of Table [J^ 



In investigating models with the HSR parameters, we impose an additional prior N > 25 
as default. Basically, subscripts in names of models represent varied parameters in the 
model, and other parameters absent in the subscripts are fixed to the default values. The 



only exception is e* which is varied in a small range [0, 10 -4 ] in M^P 6 !. This is because e* 
cannot be fixed to zero as the inflaton cannot roll when e* = 0. Regarding priors on A s , 
we adopt a common prior 2.5 < ln[10 10 A s ] < 3.5 in all models. In our analysis a fiducial 



wave number is chosen to be fc* = 0.01 Mpc~W 7 . 

We assume a flat ACDM model as background cosmology. In addition to parameters 
representing the shape of the primordial power spectra (i.e. (A s , n s , r, a s ) or (A s , e#, 77*, 
£*)), following parameters are included in the analysis as primary cosmological parameters, 

(u b ,u c ,9 s ,T, A sz ), (26) 

where Uf, and u c are the density parameters of baryon and CDM, respectively, 9 S is the 
angular scale of the acoustic horizon [62], r is the optical depth of reionization, and Asz 
is the amplitude of template Sunyaev-Zel'dovich power spectrum Cf z [63]. The priors on 
these primary cosmological parameters are listed in Table [2J 

We adopt CMB data from WMAP5 [631466] . as well as the observations at small an- 
gular scales including ACBAR [67], CBI [68], BOOMERANG [MLE] and QUAD [TJ. 
In addition, we adopt observational data of geometrical distances of the late Universe, 
including the Union data set of SN [76] , the measurement of BAO scales in galaxy power 
spectra and the SH ES measurement of the Hubble constant H = 74.2 ±3.6 [75]. A 



# 5 A similar division of models can also be found in [3U] • 

# 6 The upper bound for e* is chosen so that these models cover inflation scenarios where detection 
of primordial B-mode would be difficult with CMB observations in the near future. 

# 7 This is very near a fiducial wave number fc» = 0.017 Mpc -1 , which are suggested to be optimal for 
constraining the primordial power spectra from current data |61j . Choice of fiducial wave number is also 
discussed in [23] . 

# 8 We also performed the same analysis by adopting the data of halo power spectra from the catalogue 
of the SDSS Luminous Red Galaxies [73] instead of the BAO data. However, as long as combined with 
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models 


description 


M HZ 


HZ spectra 


M na 


0.8 <n s < 1.2, r = a s = 


M nsr 


0.8 <n s < 1.2, < r < 0.5, a, = 


M naaa 


0.8 <n s < 1.2, r = 0, -0.1 < a s < 0.1 


M e 


< e* < 0.1, 77* = £* = 


M v 


< e* < 10~ 4 , -0.1 < 77* < 0.1, = 


M er) 


< e* < 0.1, -0.1 < 77* < 0.1, = 



Table 1: Models adopted in the analysis. Mhz is the reference model that has HZ primor- 
dial power spectrum. Other models cover a model space for single-field slow-roll inflation. 
M„ s , M nsr , M nsCta have top-hat priors on the standard parameters of the primordial power 
spectra, and M e , M v , M er] have ones on the HSR parameters. Prior ranges are described in 
the last column. Regarding priors on A s , we adopt a common prior 2.5 < ln[10 10 v4 s ] < 3.5. 



parameters 


prior ranges 


u b 


[0.02,0.025] 


u c 


[0.08,0.14] 


9 S 


[1.02, 1.06] 


T 


[0.01,0.2] 


A S z 


[0,4] 



Table 2: Top-hat priors on cosmological parameters, other than parameters of the primor- 
dial power spectra. 



default data set that we conduct our analysis with is a combination of all the data above, 
which we denote as 'ALL'. On the other hand, we also make use of several other data 
sets, in order to assure ourselves that the results are not dragged by some extreme data. 
Adopted are data sets of WMAP5 alone and the combinations of WMAP5 with either 
the other CMB data at small angular scales, BAO, SN, or the H measurement, which we 
denote as +CMB, +BAO, +SN, +H , respectively. 

Computation of Bayes evidence is implemented with MultiNest [79], which is inte- 
grated in the Markov chain Monte Carlo (MCMC) sampling code CosmoMC [80] but uses the 
nested sampling [81] in stead of the MCMC sampling (See also Ref. [82]). Given a model 

data of CMB, SN and the Hq measurement, the results are very similar to those from the BAO data, 
which is also consistently observed in Ref. [74]. We suppose that this is due to the marginalization over 
scale-dependent galaxy biases. Although the matter power spectrum itself is a promising probe for the 
primordial power spectrum at small scales (See also |75j for a recent result), there are some difficulties 
in modeling nonlinear evolution of matter perturbations, biases of tracers, or effects of baryonic physics. 
For the time being, we do not adopt data of the matter power spectrum. 



S 



M, MultiNest provides chains of samples from the posterior distribution P(0|data, M) 
and Bayes evidence E(M). Convergence is diagnosed by applying the Gelman- Rubin test 
with two independent chains. Typically R — 1 < 0.01 is achieved. 



4.1 Standard parameterization 

In Table [3] we list Bayes factors for models with standard parameterization of the primor- 
dial power spectra. M ng , M ngT and M ng0ls against a reference model Mhz- What we may 
first see in the table is that almost all data sets give positive Bayes factors. Although 
there are a few exceptions given by the +Hq data set, they are by no means significant. 
This result shows us that current data negatively support a model with the HZ spectrum. 
In particular, the strongest negative support is brought from the ALL data set. Although 
Bayes factors depend on models we compare with M H z, all models we adopt, M ng , M nsr 
and M naaa , give Bayes factors larger than 2.5. On Jeffreys' scale, these results correspond 
to strong negative evidences for Mhz- Such negative preference for the HZ power spectra 
was reported by many previous studies [36 |[38Hl0l H2 | l^] and our result is consistent with 
the other recent result from WMAP5 jH] . However, significance in our result is consider- 
ably large as very recent observations are adopted. In Figure HJ we plotted the Id and 2d 
marginalized posterior distributions of the parameters n s , r, and a s from the ALL data 
set. Clearly, the negative support for Mhz originates from a poor fit of the HZ spectrum to 
the data, as the posterior distributions of n s for models M ns , M nsr and M nsOCs consistently 
show preference for n s ^ 1. 

The Bayes evidences in Table [3] also show that the current data can be well-modeled by 
scalar perturbation with power-law power spectrum without the tensor mode (i.e. M n J. 
Let us regard M ng as a base model and ALL as a base data set. When we include the 
tensor perturbation by varying r > 0, the Bayes factors decrease from 4.3 for M„ s to 
2.6 for M ngT . This suggests that current data negatively support presence of the tensor 
perturbation. In fact, the posterior distribution for r peaks at r = 0, as seen in Figured] 
On the other hand, when we include the running of the scalar spectral index, the Bayes 
factor is almost unchanged from 4.3 for M ng to 4.4 for M na0l3 . This shows that even though 
the posterior distribution of a s in Figure [T] peaks away from a s = 0, the likelihood does 
not improve much from M ng in most region of ot s ^ for the M nsCls model. We conclude 
that presence of a s is not disfavored by the data nor required in modeling the data. 

From Table [3] the preference for M ng over M ngT is commonly observed with all the 
data sets. Also M ng and M ngCls are supported almost equally since the Bayesian evidences 
for these models do not differ by 1 at most. Totally, as M ng is nested by M ngT or M ngag , 
we can safely conclude that M Ug best describes current data as far as we compare models 
with standard parameterization listed in Table [31 

The constraints on cosmological parameters for these models are summarized in Ta- 
ble HI As seen in the table, the current data can optimally be modeled by M Ug . Therefore, 



we propose the constraint for M Ug as optimal onq^j. An optimal constraint on the pri- 



^ 9 We do not derive an optimal constraint by averaging over models, as our division of models, M n 
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models 




M nar 




WMAP5 


+0.8 ±0.2 


-0.1 + 0.2 


+0.3 + 0.2 


+CMB 


+1.4 ±0.2 


+0.5 + 0.2 


+1.2 + 0.2 


+BAO 


+1.7 + 0.2 


+0.4 + 0.2 


+1.5 + 0.2 


+SN 


+1.6 + 0.2 


+0.3 + 0.2 


+1.9 + 0.2 


+H 


+0.5 + 0.2 


-0.2 + 0.2 


-0.4 + 0.2 


ALL 


+4.3 + 0.2 


+2.6 + 0.2 


+4.4 + 0.2 



Table 3: Bayes factors for models with the standard parameters of the primordial power 
spectra against the reference model Mhz- 



mordial power spectra from the ALL data set is given by 

ln[10 10 A s ] = 3.1371S (27) 
n s = 0.957+Hn, (28) 

where errors are given at 68% CL. Regarding the tensor perturbation and the running of 
the scalar index, we find no evidence at present. 

As a final remark of this section, let us comment on the choice of prior distributions 
and its relations to the resultant Bayes evidences. As we have mentioned in Section [31 
Bayes evidences depend on prior distributions. While a prior 0.8 < n s < 1.2 is often 
adopted in other studies [39|l4Tj. priors r < 0.5 and \a s \ < 0.1 may look moderately tight. 
More relaxed priors on r and a s may slightly decrease the Bayes evidences for M ngr and 
M nsCea , which however hardly affects our conclusion. 

4.2 HSR parameterization 

In Tables |5] we present the Bayes factors for models with the HSR parameters, M e , M n and 
M ev , against Mhz- From the table, first we reconfirm that the Mhz model is negatively 
supported from the data. With the ALL data set, the Bayes factors are 5.8, 4.7 and 4.2 
for M e , M v and M er) , respectively. On Jeffreys' scale, these correspond to negative support 
for the model with the HZ spectrum at decisive or at least strong level. This shows that 
models of primordial power spectra with inflationary priors can quite well explain the 
data, while the data are very poorly fitted by the HZ spectrum. 

Let us take a close look at the results so as to find plausible models of single-field 
slow-roll inflation from current data. In Figure 121 we presented the marginalized Id and 
2d posterior distributions of the HSR parameters e* and 77*. There we explicitly show the 
parameter region in e*-?7* plane that are removed by a prior A > 25 for = 0. In Figure 
[31 we also present the marginalized posterior distributions of the amplitude of curvature 
power spectrum A s and the derived standard parameters, n s , r and a s , given as functions 

M nsr and M„ sQa , is rather artificial, and not strongly based on theoretical motivations. 
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Figure 1: Id and 2d marginalized posterior distributions of the standard parameters, 
log[10 10 74 s ], n s , r and a s from the ALL data set. Models shown are M na (dashed green 
line), M nar (dotted blue line), M ns0la (dot-dashed magenta line) and the reference M H z 
(solid red line). 
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Table 4: Constraints for models with the standard parameters of the primordial power 
spectra, M Ua , M ngr and M nsQ , s , from the ALL data set. For bounded parameters, shown 
are mean values and 68% credible intervals. For unbounded parameters, 95% credible 
intervals are only shown with parentheses. 



of the HSR parameters (See Eqs. (I16H18P ). From Table [5] we consistently observe that the 
Bayes evidences for M e are larger than those for M v and M eri , which are nearly equal, for 
all the data sets we adopted. With the ALL data set, the Bayes factor of M e against M n 
and M er] are 1.1 and 1.6. Thus the current data show slight preference for M e over M v 
and M tr] , which would be significant on Jeffreys' scale. 

However, the preference for M e looks somewhat dependent on several assumptions we 
adopted in the analysis and the origin is not clear. One possible origin would be the prior 
N > 25, for it eliminates the region e* > 0.024 if rj* = 0, as can be seen from the shaded 
region in Figure HJ Thus in the models of M € , the original top-hat prior < < 0.1 is 
effectively substituted with a top-hat prior < e* < 0.024 due to the prior N > 25. As 
can be seen in the marginal posterior distribution of e* for M e in Figure |2l the effective 
prior region appears to give high likelihood, which results in boosting up the Bayesian 
evidence for M e . 

Let us examine the origin of preference for M e in more detail. We have repeated Monte 
Carlo analyses using the ALL data set, modifying the default setup in several different 
ways as follows: 

(1) A prior N > 25 is removed. 

(2) Contribution from the tensor perturbation is omitted in the power spectra of CMB 
anisotropies. 

(3) The prior ranges for e* and 77* are extended from e*, |r/*| < 0.1 to e*, |^*| < 0.2. 

(4) The higher order HSR parameter is included with a top-hat prior —0.01 < < 
0.01. 
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The results are summarized in Table [6j From the modification (1), we first note that the 
Bayes fractor for M e significantly decreases from 5.8 to 4.7 by removal of a prior N > 25. 
The resultant Bayes evidences are very comparable between models M € and M v , while 
more complicated M er} gives a smaller Bayes evidence. 

On the other hand, the largest difference between a model with large e* (e.g. M e ) and 
one with small e* (e.g. M v ) would be whether primordial tensor perturbation is generated 
or not. Therefore the preference for M e can also possibly be induced from observably large 
contribution of the tensor perturbation in CMB anisotropies. However, this possibility is 
not supported by the result from the modification (2). The Bayes factor for M e does not 
decrease by omission of the tensor contribution, or even comes to increase slightly from 6.0 
to 6.3. This is also as expected from the discussions in Section |4TT| where the model M nsT 
with nonzero r is less supported, compared with M ng with vanishing r. From Figure [3] 
we also note that the running of scalar spectral index for these models is too small to be 
observed with current data (See also constraints on a s for M nsQs in Figured]). Therefore 
the preference for M e is not from the running of the scalar spectral index either. 

With all above demonstrations, we conclude that the preference for M t originates 
from a theoretical prior N > 25. Current data is not very constraining enough for us to 
distinguish models adopted in our analysis and future observations should be await. 

We also examine the dependence of Bayes evidence on top-hat priors on e* and 77*. 
The modification (3) in Table [6] shows that doubling of prior ranges from e#, \rj*\ < 0.1 to 
e*, I?]* I < 0.2 decreases the Bayes factors for M v , M ev by 0.8 and 1.4, respectively. These 
decreases are almost as expected from the dependence of Bayes evidences on prior ranges, 
that are about In 2 ~ 0.7 and 2 In 2 ~ 1.4 for M v and M tJ] , respectively. Contrastively, the 
Bayes evidence for M e is unchanged by the modification (3) because the priors N > 25 
dominates over the original top-hat prior on e*. 

So far we restricted ourselves to models with the HSR parameters up to e and 77. Let 
us examine possible preference for models with higher order HSR parameters. We have 
included the third lowest order HSR parameter with a top-hat prior —0.01 < < 0.01 
in the modification (4). From Table [6] we observe that the Bayes evidences almost equal 
or decrease from those with vanishing £*. Therefore current data show no preference for 
models with nonzero £*, and probably models with further higher order HSR parameters, 
as well. Here we should remark that above conclusion for models with nonzero relies 
on a prior N > 25. As explicitly shown in [23], positive ~ O(0. 01) may be favored 
when we omit the prior N > 25. As we have discussed in Section [2] even for such large 
£*, a sufficient e-folding number N > 25 can possibly be achieved by higher order HSR 
parameters or other subsequent inflations. However, such complication may make the 
models less appealing. 

The constraints on cosmological parameters for the models with the HSR parameters 
are presented in Tabled From the table, we see that constraints on the amplitude and the 
spectral index of the primordial curvature power spectrum, A s and n s , are almost inde- 
pendent of priors on the HSR parameters (See also Figure [3]). Moreover, these constraints 
are identical with those of Eqs. (I27H28P from the M ng models. Therefore we conclude that 
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models M e M v M er] 

WMAP5 +2.0 ±0.2 +0.9 ±0.2 +0.6 ± 0.2 

+CMB +2.9 + 0.2 +1.7 + 0.2 +1.1 + 0.2 

+BAQ +3.2 + 0.2 +1.9 + 0.2 +1.6 + 0.2 

+SN +3.0 + 0.2 +1.9 + 0.2 +1.7 + 0.2 

+H +1.8 + 0.2 +0.6 + 0.2 +0.3 + 0.2 

ALL +5.8 + 0.2 +4.7 + 0.2 +4.2 + 0.2 

Table 5: Bayes factors for models with the HSR parameters against the reference model 
M HZ - 



models M e M v M tr} 

ALL +5.8 + 0.2 +4.7 + 0.2 +4.2 + 0.2 

(1) +4.7 + 0.2 +4.6 + 0.2 +3.7 + 0.2 

(2) +5.9 + 0.2 +4.6 + 0.2 +4.2 + 0.2 

(3) +5.8 + 0.2 +3.9 + 0.2 +2.8 + 0.2 

(4) +4.3 + 0.2 +5.1 + 0.2 +3.6 + 0.2 

Table 6: Bayes factors for models with the HSR prameters against the reference model 
Mhz- Several extensions from the default setting are investigated: (1) removal of a prior on 
e-folding number; (2) omission of the tensor perturbation; (3) imposition of extended top- 
hat priors, e*, \r]J\ < 0.2; (4) inclusion of the third lowest order HSR parameter —0.01 < 
< 0.01. 



the optimal constraint we have proposed in Section 14.11 is robust with little dependence 
on priors from different inflationary models. However, constraints on the tensor pertur- 
bations and the running of the spectral index are largely dependent on inflationary priors 
and their presence is not suggested from current data. 

5 Summary and future outlook 

We have investigated constraints on the primordial power spectra and comparison of single- 
field slow-roll inflation models using data from recent observations of CMB combined with 
measurements of BAO, SN and H . By employing Bayesian model selection, we found 
that a model with the scale-invariant HZ spectrum is strongly disfavored from current 
data, in comparison with several simple models allowing scale-dependence of the power 
spectra. We have also proposed an optimal constraint on the primordial power spectra of 
Eqs. (ETJEED. 

Adopting our somewhat artificial division of models for single-field slow-roll inflation, 
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Figure 2: Id and 2d marginalized posterior distributions for the HSR parameters, e* and 
77*, and the amplitude of the curvature power spectrum ln[10 10 A s ]. Models shown are 
M e (dashed green line), M v (dotted blue line), M trj (dot-dashed magenta line) and the 
reference Mhz (solid red line). The gray shaded region is excluded by a prior N > 25 
when = 0. 
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Figure 3: Same figure as in Figure [2] but the posterior distributions for the standard 
parameters, ln[10 10 y4 s ], n s , r and a s , are shown. 
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Table 7: Constraints for models with the inflationary HSR parameters, M ( 
from the ALL data set. For bounded parameters, shown are mean values and 68% cred- 
ible intervals. For unbounded parameters, 95% credible intervals are only shown with 
parentheses. 



M e , M v and M e7? , with a theoretical prior iV > 25, we have found the Bayes evidences for 
these models from current data are almost comparable, but there is a slight but consistent 
preference for M e over the others from various data sets. We have demonstrated several 
extension in order to identify the origin of the preference for M e . A presence of the tensor 
perturbation is not signified and we have found that the preference originates from a prior 
N > 25 which appears to make a prior ranges for e* center around the region of high 
likelihood. We have also found that higher order slow-roll parameters are not required 
from current data. Thus we conclude that while simple models of single-field slow-roll 
inflation can describe current cosmological observations, data is not enough constraining 
to distinguish these models we adopted. 

Planck [83] and several other ground-based CMB surveys, such as QUIET [84J, Po- 
larBeaR (85], etc. are now under way. In the near future, the tensor perturbation will 
be probed down to r ~ 0.01 by these surveys. This allow us to discriminate models of 
single-field inflation through the prediction on the primordial power spectra or the slow- 
roll parameters. Indeed, one of our division of models, M e , will be strongly refused, if the 
primordial tensor perturbation is not detected by those surveys. Given data from these 
future surveys, employment of Bayesian model selection would be of great help in several 
ways, such as constraining cosmological parameters in an optimal way, and assessing a 
support for a cosmological model statistically. 
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